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Abstract 

We propose a simple model of cytochrome c oxidase, including four redox centers and four pro- 
tonable sites, to study the time evolution of electrostatically coupled electron and proton transfers 
initiated by the injection of a single electron into the enzyme. We derive a system of master equa- 
tions for electron and proton state probabilities and show that an efficient pumping of protons 
across the membrane can be obtained for a reasonable set of parameters. All four experimentally 
observed kinetic phases appear naturally from our model. We also calculate the dependence of the 
pumping efficiency on the transmembrane voltage at different temperatures and discuss a possible 
mechanism of the redox-driven proton translocation. 
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I. INTRODUCTION 



The last enzyme of the respiratory chain of animal cells and bacteria, cytochrome c 
oxidase (CcO), operates as an efficient nanoscale machine converting electron energy into a 
transmembrane proton electrochemical gradient . The ATP (adenosine 

triphosphate) synthase enzyme uses this energy to synthesize ATP molecules serving as the 
"energy currency" of the cell. The process of energy conversion starts when a molecular 
shuttle, cytochrome c, delivers, one by one, high-energy electrons to a dinuclear copper 
center, Cu^, located near a positive side (P— side) of the inner mitochondrial membrane (see 
Fig. 1). In recent time-resolved optical and electrometric studies jsl of the CcO transition 
from the oxidized (O) state to the one-electron reduced form (E), a single electron is donated 
to the Cua redox center by a laser-activated molecule of ruthenium bispyridyl (RubiPy). 
Thereafter, in a few microseconds (~10 /is), a major part of an electron density (~70%) is 
transferred from the Cu^ center to the low-spin heme a (Fe-a). Heme a is located within the 
membrane domain at a distance about 2/3 of the membrane width, W, counting from the 



N-side 



lOj. Within a time interval of approximately 150 /is, about 60% of the electron 



population is transferred from heme a to heme (Fe-a^). Heme 03, jointly with the next 
electron acceptor in line, a copper ion Cu^, form a binuclear center (BNC, sites R and B 
in Fig. 1), serving as an active catalytic site for dioxygen reduction to water. The redox 
centers a, 03, and Cu^, are located approximately at the same distance (2/3 W) from the 
N-side of the membrane as heme a. 

The next phase (with a time scale of the order of 800 /is) is characterized by a complete 
electron transfer to the copper ion Cu^. Time- resolved measurements 8| show that the first 
"10 fis" phase of the electron transfer process is not accompanied by a proton transfer, but 
the "slowness" of the second "150 /is" and the third "800 /is" phases hints to the proton 
participation during phases. 

The proton path from the negative side of the membrane (N-side) toward the P-side (for 
pumped protons) and toward the binuclear center (for substrate or "chemical" protons) goes 
through the residue £'278 (for the Paracoccus denitrificans enzyme jsj). These residues 
are located at the end of the so-called D-pathway (Fig. 1). A fraction of the substrate 
protons can also be delivered to the BNC via an additional K-pathway, which we will not 
consider here. The proton to be pumped is supposed to move from E278 (schematically 
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shown as the site "D" in Fig 1) to an unknown protonable "pumping" site X (hkely a heme 



and, thereafter, via an additional protonable 



as propionate], located above the BNC 

site C 12, ll3|, to the P-side of the membrane. After a fast reprotonation from the D- 
channel, the residue £'278 can donate a substrate proton to the catalytic site near the BNC 
(probably, to an OH" ligand of Cu^ 9]). It is assumed j^, 9] that during the second 
"150 /is" phase, the first (pre-pumped) proton moves from the residue E278 (the site "D" 
in Fig. 1) to the pump site X, whereas in the third "800 yus" phase, the second (substrate or 
chemical) proton populates a catalytic site Z near the BNC. In the final phase, which occurs 
in 2.6 ms, the first proton (in X) is translocated (via C) to the P-side, which is characterized 
by a higher electrochemical potential than the N-side of the membrane. 

As a result of all these processes, two protons are taken from the N-side of the membrane, 
and one electron is taken from the P-side, and eventually one proton is pumped to the P-side. 
Moreover, one proton and one electron are consumed at the catalytic site to finally produce 
a water molecule around the BNC. It should be noted that kinetic phases with similar time 



scales (10 /is 100 fis 1000 fis) have been reveal 



transition steps between the states of the enzyme 12 



Kinetic data obtained in experiments 



12 



14 



ed in other CcO enzymes at various 



14, 



15|. 



15| reflect important details of the still 



elusive proton pumping mechanism in cytochrome c oxidase. To extract these details and 
gain a deeper insight into the operating principles of the CcO proton pum p, i t is necessary 
to compare results of experiments with theoretical predictions. In Ref. [l6|], a simplified 
empirical valence bond (EVB) effective potential was combined with a modified Marcus 
equation to model time- dependent electron and proton transfers in CcO in the range of 
milliseconds. However, this approach was applied to the single transfer event, not to the 
sequence of events, and the obtained time scale (one microsecond) differs by orders of mag- 
nitude from the experimental data (about 100 microseconds) . A computational analysis 
of the CcO energetics was presented in Refs. [17|, 



18 



m, 20, 



- 2l| with molecular models 
reproducing energetic barriers for the proton transfer steps 17, ISJ. The obtained energetic 
map of the proton and electron pathways in the CcO enzyme can be converted into a set of 
rate constants, which qualitatively explains the kinetics and unidirectionality of the pump- 
ing process. However, these studies do not result in a quantitative model of the efficient 



CcO proton pump. Moreover, the error range of these semi-microscopic calcu 



kcal/mol) is sometimes higher than the difference between the energy barriers 18l. l20|. 



ations 
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Kinetic models of the proton pumping process were also discussed in Ref. [22|. Within 
the master equation approach, it was shown that the proton pumping effect can be achieved 
in a simplified system having one redox and two proton sites and, with a higher efficiency, 
rj ~ 0.9, for the design with two redox and two protonable sites, which are electrostatically 
coupled to each other. However, this work does not contain any predictions for the kinetics 
of the pumping process in more reahstic set-ups, with at least four redox sites (Cu^i, heme 
a, heme 03, and Cu^) and two protonable sites (a residue £^278 and a pump site X). To 

i resort to a Monte 



find proper parameters for the proton pump, the authors of Ref. 
Carlo search in a multidimensional parameter space. It is hard to imagine, however, that 
a random search can provide a reasonable set of parameters which will comply with all 
physical restrictions of real pumps. In general, for comprehensive theoretical studies, it is 
preferable to determine the relevant parameters of the system using detailed microscopic 



calculations (see, e.g., Refs. Q, [ITI, ll8|]). However, the huge computational complexity of 



biological structures makes such an approach extremely difficult. In our paper, we include 
reasonable estimates for the system parameters into a model describing almost simultaneous 
electron and proton transfer processes and compare the obtained kinetics to experimentally 
observed time scales and site populations of cytochrome c oxidase js]. 



The time evolution o 
data of Refs. 



0, Q, Q, 



the proton pumping process in CcO, related to the experimental 



15|, was discussed in Refs. 



19 



20 



23j. In these works, the 



kinetics of the electron-proton system is broken down into a cascade of quasi-equilibrium 
states characterized by distributions of electrons and protons over the sites, as well as by 
a set of transition rates corresponding to specific kinetic phases. It should be emphasized, 
however, that many electron and proton transfers can be separated by only a nanosecond 
time scale, and, consequently, the experimentally observed kinetic rates comprise contribu- 
tions of several almost-simultaneous individual electron and proton transfer events 



10 



23|. 



Correspondingly, an approach taking into account the kinetic inseparability of electron and 



proton transitions can be useful for understanding recent experimental 



10 



19 



findings 



20 



We note 



23( 1 were deduced 



that kinetic coefficients used in the theoretical analysis of Refs. 
from experiments without independent microscopic calculations of the heights of individual 
electron and proton barriers. 

In the present paper, we analyze electron and proton kinetics in cytochrome c oxidase 
within a simple physical model including four redox centers and four protonable sites elec- 
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trostatically coupled to each other in the presence of a dissipative environment. Using the 
master equation approach, we reproduce all four kinetic phases observed in Ref. js] for a 
reasonable set of parameters. It should be emphasized that we have performed extensive 
numerical studies for a wide range of parameters and we found that our model of proton 
pumping is quite robust to significant variations of the parameters. The specific set of pa- 
rameters presented below gives a very good agreement with the experimental data of Ref. |8| . 
We consider a single cycle of events, which starts at t = with one electron transfer to the 
Cu^ center and finishes at the moment t = ts, when the redox site Cu^ is completely re- 
duced. Notice that the injection of additional high-energy electrons is necessary to maintain 
this nonequilibrium state of the CcO enzyme. We also determine the efficiency of the pro- 
ton pumping for our model and its dependencies on the temperature and transmembrane 
voltage. 

The rest of paper is structured as follows. Our model and its parameters are presented 
in Section II. Results of numerical studies are shown in Section III and discussed in Section 
IV. Section V contains the conclusions of our work. The detailed derivation of the master 
equations and the measurable variables is presented in the Appendix. It should be noted that 
while the results of this paper are obtained in the classical regime, our approach (based on 
quantum transport theory) can be used to examine fine quantum effects and, consequently, 
the detailed derivation is worth presenting here. 



II. MODEL 



As in the real CcO enzyme 



24 



25 



26 



27 



the redox chain of the present model 



includes four centers: Cu^ (site A), heme a (Fe-a, site L), heme 03 (Fe-as, site R), and Cu^ 
(site B), as schematically shown in Fig. 1. The transport chain for protons has four sites: D 
(presumably related to the residue E278 near the end of the D-pathway), X (the pump site 
above the BNC), a protonable site C placed on the way from the X-site to the P-side of the 
membrane, and, finally, a protonable site Z located in the proximity of the BNC and related 
to the 0H~ ligand of Cu^ (see Fig. 1). The sites B and Z serve as final destinations for the 
injected electron and for the substrate proton, respectively. We assume that the electron 
can be transferred between the pairs of redox states A and L, L and R, R and B; and that 
protons can be translocated between the pairs of protonable sites D and X, X and C, as 
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well as D and Z. 



To provide an "openness" of the CcO enzyme, which is inherent in the living systems [22| , 
we allow proton transitions between the site D and the negative side of the membrane as well 
as between the site C and the positive side of the membrane. Protons are delivered to the 
catalytic site Z partially through the K-pathway . This channel can be incorporated 
into our model, but, for simplicity, it will be neglected. The N- and P- sides of the membrane 
play roles of proton reservoirs which work as a source (N-side) and a sink (P-side) of protons 
for the enzyme. The redox sites are disconnected from electron reservoirs, and only one 
electron is injected into the redox chain at the initial moment of time, t = 0. 

With the condition of single-occupation of each individual site, the system can be pop- 
ulated with up to four protons. Following the setup of Ref. j^, we assume that CcO is 
populated with a single electron initially located on site A. To quantitatively describe this 
system we introduce 64 basis states /i = 1,...,64 (see Appendix). The time evolu- 
tion of the probabihty distribution over the basis states, {pit)) = {(p^(t))}, is governed by 
the system of master equations, Eq. (1A32I) . with the solution given by Eq. flA34l) in the 
Appendix. The time-dependent probability distribution (p(t)) allows us to determine the 
average populations of all electron and proton sites, (wq,) and {Np), as functions of time. We 
can also calculate the number of protons, {Np(t)), translocated to the positive side of the 
membrane [see Appendix, Eq. flA36l) ]. The value of (A^p) taken at the end of the pumping 
cycle [t = Ib) determines the pumping efficiency t] defined j2J] as the number of protons 
pumped across the membrane per electron consumed: 

yi={N^{tB)). (1) 

Note that the efficiency rj can take negative values in the case when protons move back from 
the positive side to the negative side of the membrane. 



A. Electrostatic interaction 



The electrostatic interaction between the redox (a = A, L, i?, B) and protonable {j3 = 
D,X,C, Z) sites plays a pivotal role in the electron-proton energy exchange. It should be 
noted that we consider here only direct Coulomb interactions between electron and proton 
subsystems and between protons themselves. This removes strict geometrical restrictions 
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on the relative positions of electron and proton active sites imposed in our previous model 



291] based on the Forster-type energy exchange between electrons and protons. Microscopic 



calculations of the electrostatic parameters, u^is and M/3/3', involved in the Hamiltonian Hq [see 
Appendix, Eq. flAip ]. require a detailed knowledge of the CcO structure complemented by 
the comprehensive dielectric map of the enzyme 17|, |30|, l3l|] • Instead, we tune the Coulomb 
energies to get the best possible fitting of the time scales and site populations measured in the 
experiment Q| . The obtained values of Coulomb parameters correlate well with information 
about the distances between the active sites 



24 



25 



261, 



27 



28l | for reasonable values of the 



effective dielectric constants. 

To describe the experimentally observed kinetic phases of the pumping process, we assume 
that the coupling, ubz = 630 me V, between the copper ion Cu^ and the catalytic site Z 



(likely an OH ligand of Cu^ 



]) and the coupling, urx = 555 meV, between heme 



03 and the pump site X are higher than the electrostatic energies urz = 530 meV and 



24 



25|, 



26 



27 



28l | performed 



ubx = uxz = 510 meV. Structural studies of the CcO enzyme 

at a resolution of about 2 A show that the BNC redox sites R (heme 03), B (Cub) and the 
protonable sites X and Z are separated by a distance of the order of 6 A. The value of the 
electrostatic coupling between these sites, u ~ 600 meV, roughly corresponds to the effective 
dielectric constant, e ~ 4, which is of frequent use for a description of a dry protein interior 







20 



211 ]. It should be emphasized, however, that the concept of dielectric constant is 



28 



site D) and the sites L and R are 



321 ]. We estimate the electrostatic 



not completely appropriate for a calculation of Coulomb potentials in the heterogeneous 
environment inside and near the BNC jsol . 31]. 

The distances, rLC),rRD, between the residue E27S 
almost the same: = 12.3 A, r^/) = 12.8 A 
coupling between these sites as uld ~ urd — 75 meV, which corresponds to the higher 
dielectric constant e ~ 15. We consider a smaller dielectric constant, e ~ 10, for the inter- 
action, ulx = 100 meV, between the sites L and X separated by the distance r^x ^ 14.2 
A 2o| . Distant-dependent dielectric constants, e(rQ^), are common in protein electrostatics 



17, 



30 



3l|. 



Note that here, as in the models of Refs. [o], Q, the electrostatic coupling, urx, between 
heme 03 (site R) and the site X is stronger than the interaction, u^x, between heme a (site 
L) and the pump site X. For the other parameters we choose the following values (in meV): 
udx ~ 60, udz ^ ubd ^ uxc ~ 100, uad ~ 25, uaz ~ 20. The Coulomb energies 
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ucz, udci uax, uaci urc, ulci ci-nd ulz are assumed to be near 30 meV. Despite the fact that 
these energies are about or higher than the temperature energy scale, T = 298 K ~ 26 meV, 
they have a minor influence on the performance of the model. 



B. Energy levels of the sites 

We assume that the difference A/i ( ]A37[) between the electrochemical potential /ip of the 
P-side and the potential /in of the N-side of the membrane is about 210 meV at standard 
temperature, T = 298 K, with yUp = 105 meV and /^n = —105 meV. This corresponds to 
voltage V ~ 150 meV applied across the membrane. We include the electron charge in the 
parameter V and measure voltage, along with other energies, in units of meV. According to 
Eqs. (IA38I) . the energy levels, and ep, of the electron and proton centers are shifted from 
their intrinsic values e^a^ and e^^^ depending on the voltage V and on the positions XaiXjj 
of the active sites. To estimate the electron and proton energies, we take into account the 



facts [g] that cytochrome c delivering electrons to the CcO enzyme has a redox potential 
of order of 250 meV, and that the total drop of electron energy between cytochrome c and 
the dioxygen reduction site B is about 550 meV. The equilibrium midpoint potentials 

m 

of the Cua center {Em — 250 meV) and heme a {E^ — 270 meV) can also be used as a 
general guide for estimating energies |33], although the real parameters can deviate from 
the estimated values. 

We find that our model performs with the high efficiency, rj ~ 0.95, and reproduces all 

I I 

experimentally observed kinetic phases |8i] for the following set of electron intrinsic energies 
(inmeV): e^^ = —175, ef^ = —240, = —185, e^-* = —155, and for the following energies 
of protonable sites (in meV): e^"* = —100, e^^ = 250, e^?-* = 195, and e^^^ = —65. It should 
be noted that in the presence of the transmembrane voltage, \^ = 150 meV, the electron 
energy levels of A and L sites, Ea = — 250, = —265, are close to the values extracted from 



equilibrium redox titrations (see also Ref. 3J], where an estimation, {ea — Sl) — 18 meV, has 



been obtained). For energies of other redox sites we use the values: Er = —210, Eb = —180. 
The energies of the protonable sites are also shifted with voltage, V = 150 meV, present: 
Ed = —85, Ex = 295, Eq = 270, and Ez = —40. It should be stressed that the energy, ex, 
of the pump site X is set to be higher than the potentials of the proton reservoirs on both 
sides of the membrane: ex > yUp > yUN- However, the presence of an electron on the site 
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R decreases the proton energy to the level, e ~ {sx — urx) ~ —260 meV, which is below 
the energy of the D-site and below the electrochemical potential, /xn = —105 meV, of the 
N-side of the membrane. As a result, the pump site X is populated with a pre-pumped 
proton. When the chemical proton moves to the site Z and the electron is transferred to the 
-B-site, the energy level of the X-site returns to the initial position, Sx = 295 meV, since 
the electron and proton charges of the catalytic site compensate each other, ubx = uxz- 
The high-energy pre-pumped proton can now move to the site C and, after that, to the 
P-side of the membrane characterized by the potential /xp = 105 meV. A large energy gap, 
{sx — £d) ^ 380 meV, significantly suppresses the return of the X-proton to the site D and 
to the N-side of the membrane. 



C. Reorganization energies and transition rates 



Part of the energy delivered to the redox center Cua at the initial time, t = 0, is dissipated 
to an environment characterized by sets of electron {Xaa') and proton (A/3/3') reorganization 
energies. To be efficient, the proton pumping process should occur with minimal energy 



dissipation. It is shown in Ref. 



35I ] that the reorganization energy for the a to electron 



transfer in the CcO enzyme can be as low as 100 meV. Similar estimates apply for the proton 
reorganization energies {sgI. Is?!. Here, we use the higher energy parameter, Xal = 200 meV, 



for the A-to-L transfer and accept the lower value, A 



and proton transitions. It is argued in Refs. 



34 



38 



A 



/3/3' 



100 meV, for other electron 



39| , that for the Cua heme a electron 



transition the reorganization energy must be in the range from 150 meV to 500 meV, which 
is much lower than the typical values of the reorganization energy for electron transfers in 
protein. The low values of electron reorganization energies (A ~ 2 - 4 kcal/mol) have also 
been calculated for electron transfer reactions in Rhodobacter sphaeroides 30[ |. 

To reproduce the initial kinetic phases, we use the following tunneling energies: Aal ~ 
0.9 /ieV, and Alr ~ A^i?, ~ 14 yueV. The parameters Alr and Arr describe the electron 
transfers, which are coupled to the slower proton transitions characterized by the energy 
scales: Adx ~ Acx ~ 0.3 yueV, and Azd ~ 0.2 fieV. It should be noted that the electron 
transfer between heme a and heme 03 can occur in a nanosecond time scale 401]. The 
hydrogen-bonded chains in proteins are also able to conduct protons in nanoseconds or 



faster 



41 



421. 
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We also select the values Fn ~ Tp ~ 17 ms ^ for the parameters Fn and Fp, which 
determine the flow of protons through the enzyme. These parameters Fn and Fp are of the 



same order as some of the transition rates used in Ref. 



221. 



III. RESULTS 



A. Four kinetic phases 

In Fig. 2, starting at t = 0.1 /is, we show a process of population and depopulation of the 
electron, nA,L,R,B, and proton, Nd,x,c,z, sites as well as the time dependence of the average 
number of protons pumped to the positive side of the membrane, Np. From here on we drop 
the brackets (...) denoting the averaging over the environmental fluctuations and over the 
states of the proton reservoirs. The calculations are performed for the standard conditions 
[fip = 105 meV, /in = —105 meV, ApH = — 1, T = 298 K) and for the transmembrane 
voltage V = 150 meV. We assume that initially a single electron is located at the site A 
(Cu^), and a proton occupies the site D. This means that at t = only one element of the 
density matrix is not equal to zero: p2(0) = 1. 

During the first phase of the process the electron moves from the site A to the site L 
(heme a). In ~ 10 /is near 70% of the electron density is transferred to the heme a (site L) 
with the remaining 30 percent distributed almost equally between the site A (Cua) and the 
site R (heme 03). This corresponds roughly to the 70 percent electron population of heme 
a after the first 10 /is phase observed experimentally in Ref. j^. No pronounced changes in 
populations of the protonable sites accompany this stage [see Fig. 2 (b)]. 

The second phase of the electron transfer is postponed by the time 150 /is, despite the 
fast intrinsic transition rate between the L and R redox sites. Besides the 55 meV potential 
difference between the sites R and L, the electron transfer in this phase is hampered by the 
involvement of protons. It is evident from Figs. 2a and 2b that, with a microsecond delay, 
the slightly uphill electron transfer from the site L to the site R is followed by the proton 
translocation from the site D {Ed = —85 meV) to the pump site X having much higher 
initial energy. Ex = 295 meV. This transition has been made possible by the strong R-X 
Coulomb attraction {ujix = 555 meV) lowering the effective energies of both electron and 
proton sites. In line with the experimental data ^] at t = 150 /xs, the electron density is 
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located mainly on the site R (60%) and partially on the sites L (~20%), and on the site 
B (~15%). The site A is practically empty at this stage. It is important that at almost 
the same moment of time {t ^ 150 /is) the population of the protonable pump site X also 
reaches its maximum (~65%). 

It is evident from Fig. 2b that the occupation of the pump site X is accompanied by 
the monotonic population of the the protonable catalytic site Z, thus lowering the energy 
of the redox site B from its initial level, 63 = —180 meV, to the final value of the order of 
—820 meV (see also Fig. 3). The population of the 5-site, n^, closely follows (with a small 
delay) the population Nz of the proton catalytic site Z (see Figs. 2a and 2b). It can be seen 
from Fig. 2b that in ~300 microseconds the pumped proton moves from the site X to the 
transient site C, placed between X and the P-side of the membrane, and after that to the 
positive side of the membrane. 

In the third phase (t ~ 1 ms), the substrate (chemical) proton (Fig. 2b) occupies the 
catalytic site Z, Nz > 0.8. Then, with a microsecond delay, the electron (Fig. 2a) is 
transferred, ub > 0.8, to the i?-center (Cu^), so that the heme a is practically re-oxidized, 
ul ~ 0.02. This stage is correlated with the 800 /is phase mentioned in Ref. 81]. 

In the fourth phase (t ~ 3 ms), the pumped proton (Fig. 2c) moves to the positive site 
of the membrane, Np ~ 0.95, the substrate proton populates the site Z, Nz = 1, and the 
electron is almost completely transferred to the site B, ub — 1. On average, about 1.3 
protons are taken from the N-side of the membrane during the whole process. 

The variations of the average electron energy, E^i = (i^ei), and the total energy of the 
system, 

Etot = (Ho) + ijpNp + /iNiVN, 

with time are shown in Fig. 3. Here Hq is the basic Hamiltonian of the system (lAip . H^i 
is the Hamiltonian of the electron component flA2|) . A^p and An are the average numbers 
of protons (1A36I) translocated to the P- or N-side of the membrane, respectively. At the 
beginning, the electron has energy 

^ei(O) ~ {SA - uad) ^ -277 meV, 

and at the end of the process its energy sinks to the level 

Eei{5 ms) ~ (eb - ubz - ubd) - -828 meV 
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with the total drop AE^ — 551 meV, corresponding to the experimental value [6j. The total 
energy of the system, _Etot, shows a decrease of the order of Ai^tot — 271 meV, which is less 
than the drop of electron energy since one proton gains the energy during its pumping to 
the positive side of the membrane. 



B. Pumping efficiency 

It follows from Fig. 2(c), that at t = = 5 ms, the average number of pumped protons, 
A^P, reaches its peak value, which can be used as a definition 2^ of the pumping efficiency rj: 
7] = Np{tB)- According to this definition, the present model demonstrates an almost-perfect 
performance with an efficiency n ~ 0.95 at T = 298 K, Au = 210 meV, V = 150 meV. 

n n 

This is comparable to the efficiency of cytochrome c oxidase [3|, 16|] pumping one proton 
across the membrane per one electron consumed at the oxygen reduction site. We find 
that the definition of the efficiency t] introduced above is not sensitive to the choice of the 
specific moment ts = 5 ms, since the number of pumped protons, Np{t), does not decrease 
noticeably with time during the interval from 5 ms to more than 100 ms at the standard 
conditions. 

In Fig. 4 we plot the pumping efficiency rj versus the transmembrane voltage V at three 
different temperatures: T = 150 K (blue dashed line), T = 298 K (green continuous curve), 
and T = 450 K (red dash-dotted line). We assume that the electrochemical gradient A/i 
varies in accordance to Eq. (]A37|) where ApH = — 1. AtT = 150K the pumping efficiency r) 
is almost constant at low voltages, V < 150 meV, with a subsequent drop at high voltages. 
The pump works better at room temperatures, T = 298 K, and keeps the efficiency steady up 
to voltages V ~ 200 meV. Notice that in this case the efficiency rj, which is proportional to 
the average number of pumped protons, becomes negative at > 270 meV. The performance 
of the model is significantly deteriorated at high temperatures, T = 450 K, when the proton 
flow is reversed starting with the relatively low voltage gradient, V ~ 180 meV. 



IV. DISCUSSION 



The obtained time evolution of the electron and proton populations (see Fig. 2) features 
four experimentally observed phases of the proton pumping process: the flrst "10 /is" phase. 
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when the electron is transferred from Cua (site A) to heme a (site L); the second "150 /is" 
phase when the electron moves from heme a to heme 03 (site R), and, with a microsecond 
delay, a proton partially occupies the pump site X; the third "1000 /is" phase when the 
"chemical" proton is transferred to the catalytic sites Z and, a slightly later, the electron 
is transferred to the ultimate electron acceptor Cu^. In the fourth phase, at t ~ 3 ms, the 
pre-pumped proton is released to the P-side of the membrane. 

It should be emphasized that, contrary to the models proposed in Refs. , this 

process cannot be described as a sequence of transitions between clearly defined quasi- 
equilibrium states since many electron and proton transfers occur in a very short time 
one after the other. The present theoretical model, which includes four redox sites (two 
copper centers and two hemes) and four protonable sites, is able to explain the efficient 
performance {r] ~ 0.95) of the real cytochrome c oxidase pumping almost one proton per 
one electron consumed against the electric potential difference, V > 150 meV, and against 
the transmembrane electrochemical gradient, A/i > 210 meV. We stress that all four kinetic 
phases appear naturally in our model for a reasonable set of the system parameters without 
artificial inclusions of consequent transfer processes. 

The mechanism of the proton pumping analyzed above is based on the direct electrostatic 
interaction between the redox and protonable sites, especially between the electron located 
on the site R (heme 03) and the proton located on the pump site X. The Coulomb coupling 
between the redox site Cub and the protonable catalytic site Z plays a very important role 
as well. The proton to be pumped is sequentially translocated to the P-side of the membrane 
from the sites X and C. At the beginning of the process these sites are empty since their 
energy levels are assumed to be higher than the energy levels of the proton source (/iN and 
Ed) and the proton drain (yUp). After the first "10 /is" phase the energy level of the L-site is 
slightly (~ 55 meV) lower than the energy level of the i?-site. However, an interaction with 
the environment facilitates the slow electron transfer to the site R. The population of the 
site R with the electron is accompanied by the lowering of the X-site energy level followed 
by the proton translocation from the site D to the pump site X. Because of the strong X-R 
electrostatic attraction, the effective energy of the i?-electron drops below the energy of the 
L-site, which results in the second "150 /is" phase where the major part (~ 60%) of the 
electron density is concentrated on the site R, and the pump site X is partially (~ 65%) 
populated with a proton. The electron transfer to the site R also leads to lowering the 
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energy of Z-site, thus inducing a monotonous population of the catalytic protonable site Z. 
No switch redirecting protons to the site D or to the site X (as proposed in Ref. 43|]) is 
needed here because both of these sites can be populated from the site D. 

It should also be emphasized that these three processes: the electron transfer to the R- 
site, the occupation of the pump site X, and the translocation of a proton to the Z-site, are 
strongly correlated in time. The proton transfer to the Z-site digs a deep potential well for 
the electron at the site 5, and in the third (~ 1 ms) phase the electron falls into this well. 
Afterwards, the Coulomb attraction between the pre-pumped X-proton and the electron 
is almost compensated by the electrostatic repulsion between X and Z protons, and the 
energy level of the X-proton returns to its original high value. The reverse translocation of 
the X-proton to the site D is strictly suppressed since now the energy difference between 
the sites X and D {Ex — Eo ~ 380 meV) significantly exceeds the reorganization energy 
Xox as well as the temperature broadening, l^fX^x^ ■, of the transition rates in Eq. (]A27p . 
However, the pre-pumped proton can easily move to the slightly (~ 25 meV) lower energy 
level Eq^ and, after this, to the positive side of the membrane characterized by the even 
lower electrochemical potential /ip = 105 meV. Our model does not require any nonlinear 



gates 



431] to prevent a proton leakage from the positive to the negative side of the membrane. 



V. CONCLUSION 



We have analyzed a simple model describing the kinetics of the proton pumping process in 
cytochrome c oxidase initiated by a single-electron injection. Within our model, this electron 
is subsequently transferred along four sites electrostatically coupled to four protonable sites. 
We have shown that the energy loss by this electron facilitates the proton transfer against 
the transmembrane voltage from the negative to the positive sides of the membrane with 
the efficiency ?7 ~ 0.95. In contrast to previous studies, we have not broken the electron 
and proton transfers into a series of transitions between the independent quasi-equilibrium 
states but examined inseparable dynamics of the pumping process. We have derived the 
master equations of motion and solved them numerically for a reasonable set of the system 
parameters. The obtained time evolution naturally encompasses all four experimentally 
observed kinetic phases. 
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APPENDIX A: MASTER EQUATIONS 



The kinetics of charge transfer in the CcO enzyme can be described by a set of master 



equations 
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44 



45j. For completeness we present here a derivation of these equations. 



29 
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471 ]. even though in this paper 



We start from the formahsm of second quantization 
we only discuss the classical results, with an examination of quantum coherent effects to 
be performed in the future. Electrons, located on the redox sites a [a = A, L, R, B), are 
described by the creation and annihilation Fermi operators 0^,0^, and protons located on 
the protonable sites (3 {(3 = D,X,C, Z) are described by the creation and annihilation 
Fermi operators 6/3. The spin degrees of freedom are neglected; thus, each site can only 
be occupied by a single particle. The electron population of the a-site, n^, is expressed as 
= o-t^a, and for a proton population Nf^ on site (3, we have the relation: = h'^hp. 
Protons on the negative (N) and on the positive (P) side of the membrane (cr = N,P) are 
continuously distributed over the space of an additional "quasi-wavenumber" parameter q 
and characterized by the creation and annihilation Fermi operators d^^, dqa with the density 
operator A^^^ = (i+^rfg^. 



1. Hamiltonian of the system 

The total Hamiltonian H of the electron-proton system incorporates a basic term, 

Ho = ^ EaUa + ^ s pN p + ^ uppiNpNp, - ^ Uc^pn^Nn, (Al) 

where the first and second terms describe the electron (a) and proton {(3) sites with energies 
Ea and Ep, respectively, and the third and fourth terms are responsible for the Coulomb 
interaction of protons with each other and the electron, respectively. It should be noted 
that in our single-electron model there is no inter-electron Coulomb interaction. We will 
also calculate the energy of the electron component, which is determined by the Hamiltonian 

= X] ^"'^^ ~ X] ^apNpria. (A2) 

For protons in the N-side and P-side reservoirs we introduce the Hamiltonian 

H^P = ^Eg^Ng,, (A3) 
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with the energy spectrum Sqa, whereas proton transitions between site D and the N-side of 
the membrane, and site C and the P-side are given by the transfer Hamiltonian 



characterized by the coefficients T^n and T^p. The component 



(A4) 



tun 



/3/3' 



(A5) 



is responsible for electron tunneling between the pairs {aa') of sites A-L, L-R, R-B, and 
for proton transitions between the pairs of sites D-X, X-C, and D-Z, with the cor- 

responding amplitudes A^a' (for electrons) and Aj3j3i (for protons), where A^^, = Aa'a and 



Protons are delivered from a solution to the active site D by the water-filled D-channel 



5|,l7|. It was argued [13, 



41 



48l | that the D-channel contains hydrogen-bonded chains of water 



molecules, which can convey protons via the Grotthuss mechanism. In this case, the proton 
transfer is considered as a collective motion of a positive charge through the chain, but not 



as a motion of an individual proton. According to another point of view (see Refs. 
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49 



50 



5l|). the proper orientation of water molecules required for the Grotthuss mechanism is 



characterized by a much smaller energetic penalty than the electrostatic barriers associated 
with a proton transfer through the channel. Thus, the dominant contribution to the kinetic 
rate of the proton transport in proteins is provided by the electrostatic energy. In the present 
work we model proton transitions (between the site D and the N-side as well as between the 
site C and the P-side of the membrane) by the Hamiltonian Htr (1A4I) with matrix elements 
that do not specify the transfer origin. 

The transport of protons between the active sites D-X, X-C, and D-Z are described 
by phenomenological coefficients A/j/?/ in the Hamiltonian (lASp . To obtain kinetic rates 
for proton transitions between active sites in the presence of an environment, we resort 
to the Marcus formulation of the problem. The relevant approach based on the empirical 



valence bond method 



52l | has been developed in Ref. [5l|]- As shown in Refs. [la, ll7|, l49| . 



the modified Marcus relations can be successfully applied for modelling the proton transfer 
steps in cytochrome c oxidase. 
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2. Environment 



To take into account the interaction of the electron-proton system with its environment, 
we introduce a term Hf,„- 



Hp 



(A8) 



j ^ 3 \ a p a / 

where A^o- = 'Yliq^qa is the total number of protons in the a-reservoir (a = N, P). The 
environment is represented as a set of harmonic oscillators 53|, l5J] with coordinates Xj, 
momenta pj, masses rrij, and frequencies ojj. The shifts Xja^ xjp, and xja define coupling 
strengths of electrons and protons to the environment. The total Hamiltonian H is the sum 
of all above-mentioned components: 

H = Hq + H^p + Htr + Htun + -f^env (A7) 

With the unitary transformation, 

K j \ a /3 a / 

the total Hamiltonian, H' = U^HU, can be transformed to the form 

qa j 

- ^ Tq^d^^ho - J] T*^b+dgt^ - ^ Tgpd^pbc - T*pb+dqp, (A9) 
where the operators, 

i 

Q/3/3' = <5j/3 = A/3/3' exp{i - X^^/)}, (AlO) 

i 

describe the effect of the environment on the electron and proton transitions. The protonable 
site C is located near the P-side of the membrane, and the site D is tightly coupled to the 
N-side by the D-channel. It is reasonable to assume, therefore, that C-to-P and N-to-D 
proton transitions have a negligible effect on the equilibrium position of the j-oscillator 
of the environment: Xjc = Xjp, Xjd = a^jN, so that the corresponding phase factors in 
Eq. (lAlOp related to the Hamiltonian -fftun — U^H^^^-JJ and to the total Hamiltonian (1A9P 
can be omitted. 
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3. Basis states and eigenenergies 



To quantitatively analyze the system with a single electron and with up to four protons we 
introduce a basis of 64 eigenstates of the Hamiltonian Hq: |1) = aj^|0), |2) = a\b~^\0), |3) = 



a+6+ |0), |4) = a+6+|0), |5) = a+6+|0), |6) = a+6+6+ |0), . . . , |64) = a+b+b+b+b+\0). Here |0) 



is the vacuum state of the system with no electrons and no protons, |1) is the state with an 
electron on site A, |2) is the state with an electron on site A and a proton on site D, |3) has 
one electron on site A and a proton on site X, |4) describes the state with an electron on 
site A and a proton on site C, and so on. Finally, |64) is the state with a single electron on 
site B and with one proton on each site D,X,C, and Z (i.e., a total of four protons). The 
state |1) has the eigenenergy Ei = ea, the state |2) has the energy E2 = Ea + Ed — uad, 
and the last state |64), fully loaded with four protons, has the energy 



/3=Z 

-^64 = Eb +^{ef3 - ubp) + udx + udc + udz + uxc + uxz + ucz- 

f3=D 



The Hamiltonian i^o ( 1A1I) is diagonal in the new basis: 



64 




(All) 



Other operators may have a non-diagonal form in the new basis, for example. 






(A12) 



where 



(A13) 



Here indices /i and u sweep all integers from 1 to 64. 



4. Electron and proton transitions 



In addition to the diagonal parts Hq and -ffNP, the total Hamiltonian of the system H 
contains the term responsible for the proton transitions between the N-side of the membrane 
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and the site D, and between the P-side and the site C: 

Htr = - '^iTqjsibD;fMudg^ + Tgpbc-^^d^p) p^^ + h.c. (A14) 

as well as the off-diagonal term -fftun describing the tunneling of electrons and the transfer 
of protons between the active sites, 

^tun = - 5^ A^Pm!. - 5^ Put^A^^, (A15) 

Here the operator A^v is represented by a linear combination of the bath operators 
Qal, ■■,QzD [see Eqs. (lAlOl) ]. multiplied by the non-diagonal (/i 7^ u) transition matrix 
elements (ajjai)^i,, .., (fo^fo^,)^!,: 

Afj,u = Qal{ + Qrb{ 

+ QDx{hlhx)^u + Qxc{h\hc)i,, + QDz{hlhz)^.u. (A16) 

It should also be noted that operators of the N and P proton reservoirs, dq^ and (i^p, cannot 
be completely expressed in terms of the basis operators p^y. 



5. Derivation of the master equations 

A probability (p^) to find the electron-proton system in the state |p) is determined by 
the diagonal operator p^ = |p)(p| averaged over the states of the environment and over the 
distributions of protons on both, N and P, sides of the membrane. The time evolution of 
the operator p^ is governed by the Heisenberg equation 



[p^, Htr]- - '^{A^,uP^,u - Au^^Pu^i} + ^{A^uPfMu - A„^pu^,}^- (A17) 



To derive a master equation for the probabilities (p^), we have to average Eq. ( 1A17P and 
calculate the correlation functions {A^^p^u) of the environment operators A^u flA16l) and 
the operators of the system p^^. The transition coefficients, Aq,q,/, and Tg^. are supposed 
to be much smaller than the energy scales given by the basis spectrum E^j_ (p = 1, • • • , 64). 
This means that the effective interactions with the N and P proton reservoirs [see Eq. (lA14p ] 
and with the bath of oscillators (see Eqs. flA15p . flA16p ) can be treated as a perturbation. In 
the framework of the theory of open quantum systems proposed in Ref. IsS] the correlation 



19 



function {A^^p^v) (with 7^ z/, no summation over /i and v) can be written in the form 



■/lU > 

(0)- 



+ I rfti(z[4°j(t),X'i'(ii)]-)(P.v(ii)PM.W)^(i-ii)- (A18) 



Here is a variable of the free environment (with no couphng to the electron-proton 

system), and 6{t — ti) is the Heaviside unit step function. We introduce the following 
notations for a cumulant function of two operators A{t) and B{t): 

{A{t),B{t')) = {A{t)B{t')) - {A{t)){B{t')), 

and for a commutator: 

[Ait),B{t')]^ = Ait)B{t') - B{t')A{t). 

In Eq. (lAlSp we take into account the backaction of the bath in a contrast to the approach 
of Ref. [29|] where this backaction is not included into consideration. Due to significant 
decoherence effects, off-diagonal elements of the density matrix, {p^^it)), disappear very 
fast. Because of this, the first term in the r.h.s. of Eq. flAlSI) can be neglected despite the 
non-zero value of the average unperturbed operator {A^^l{t)). The times t and ti involved 
in the integrands of Eq. flAlSp are separated by the correlation time of the correlators, 
which are similar to the function {Q Aiif) 1 Q Aiiti)) ■ The timescale Tc is determined by 



the reorganization energy \al and temperature T: Tc ~ h/ \^\alT (see Refs. [M, l56| and 
Eq. flA23l) below). We assume that transitions between the active sites have a negligible 
effect on the time evolution of the operator p^j^ between the times t and ti, which are 
separated by the correlation time Tc- Thus, the correlation functions and commutators of 
the operators p^yit) and Pu^iti) can be calculated using free-evolving functions: 

Pi,v{t) = Pi,v{h) exp{zu;^^(t - ti)}, 

where uj^j_i, = — E^j. For the correlator ( lAlSI) we obtain the formula 



{{A^S'(tlAf^\ti)){pM - {Af^\tilA^^{t))M))}- (A19) 
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With Eq. (IA16I) . we can express the cumulant {A^^}{t), A^^}^ [ti]) in terms of cumulant 
functions of the unperturbed bath operators Q^^li ■ ■ ■ i Qzd- 

= i(«>L),.r(g!Ji(t),QiT(ti)) + i(4a«),.r(QS(t),QiT(ii)) 

M{b^cbcU'{Q%{t),Q^^lt{U)) + \{bibz),A^ 

The correlation function {A^/^u~^ (ti) , A^/^u (t)) has a similar form, with cumulants 
(QALit), QAL^ih)), being replaced by (Qf^ih), Q^Aiit)), . Using the definitions 
(lAlOp of the bath operators we can calculate their correlation functions. In particular, 

{QfLit),QTiti)) = |A^,rexp{-^>V«(t-tO}exp{-Wj2(t-ti)}, 
{QT{h),Qfl{t)) = |A^,pexp{.W«(t-ti)}exp{-W(2(t-tO}, (A21) 



where 



Ml(^) = Yl ~2^^^^^ ~ ^^^^^ sin^j-r, 
i 

>VS(^) = E ^ ix,A - x,,ni - cosu;,r), (A22) 

and T is the temperature of the environment {ks = 1). These expressions can be simplified in 
the high-temperature limit when the thermal fluctuations are much faster {ujjT ^ 1) than the 
environment modes coupled to the charge transfer [5^]: W^^|(r) = A^^r, Wjf^(r) = XalTt'^, 
and, correspondingly, 

{QTiti),Qf,{t)) = |A^,|V^-(*-*^)e-^-^(*-*^)^ (A23) 
We introduce here the reorganization energy, 

A.. = 5:!^^^%^. (A24) 

j 

corresponding to the electron transition from the site A to the site L. Similar parameters 
can also be introduced for other electron transitions: from L to R, from R to B, as well as 
for proton transitions between sites D and X, X and C, and between D and the catalytic 
site Z. 
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After a sequential substitution of Eqs. (lA23p . (lA20p . (1A19P into the averaged equation 
flA17p . we obtain the contribution of the inter-site transfers into the master equation 



(A25) 



where the combined rate Hf^u contains contributions of all possible electron and proton 
transitions, 



(A26) 



The rates corresponding to the specific electron transfers, {i^al)hu, {i^lr)^iu, {i^rb)^^, and the 
rates related to the proton transfers, {KDx)fMu, {.i^xc)nv^ {.f^Dz)nv^ are all determined by the 
Marcus equations 54, l56| with coefficients given by the appropriate transition matrices. In 
particular, 

{E^ — Ey + XalY 



[i^al) 



fj,u 



TT 



XalT 



{\{a\aL)i„y\'^ + \ia\aL),y,^\'^) exp 



, (A27) 



[f^DZ, 



\Anz\' 



TT 



^DzT 



{\{bt,bz),.\' + \ibibzU\')exp 



{E^ — Ey + Xpz) 
4:XdzT 



. (A28) 



It should be noted that the ratio between the transposed rate coefficients is equal to the 
Boltzmann factor, as 



[l^AL, 



exp — 



Efj, — Ey 

T 



\I^AL)ufi 

which results in the Boltzmann distribution for the equilibrium density matrix of the system. 

The contribution, {—i[p^, i^tr]-), of proton transitions between the site D and the N-side 
of the membrane and between the exit site C and the P-side of the membrane to the master 
equation flA25P can be calculated with the methods of quantum transport theory 
The coupling to the proton reservoirs is described by the relaxation matrix. 



29|,|46|,|4 



+ Tp {|6c;^..r[l - Fpicoy^)] + \bc-u^\^Fp{iu^y)}, 



(A29) 



where the energy-independent coefficient. 



(A30) 
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determines the rate of a proton delivery to the D-site {a = N) or the rate of a proton removal 
from the C-site (a = P). We assume here that protons on the cx-side of the membrane are 
described by the Fermi distribution, 

(A31) 

characterized by a chemical potential 

As a result, we obtain the system of master equations for the probabilities (p^), as 

where the inter-site rates k^^ are determined by Eqs. (lA26[) . (lA27p . and the relaxation 
matrix, 7^,^, is given by Eq. flA29l) . 

6. Algebraic solution of the master equations 

Determination of the time-dependent solution of the master equations ( 1A32I) can be 
reduced to a purely algebraic problem. To accomplish this, we rewrite the equations flA32p 
in the form 

(^) = -EV(P-)' (A33) 

with a total relaxation matrix A^^,, where A^,y = —{n^v + liiv) at p 7^ z/, and A^^ = 
Ylivi'^i^iJ. + l^ij)- The vector (p) with the elements (p^) (p = 1, ..,64) can be represented as 
a sum of the steady-state part, p°, and the time-dependent deviation pit), as (p) = + p. 
Both the total probability vector (p) and its steady-state value satisfy the normalization 
condition: J^fiiPfj) — Yl^iP^n — ^- The steady-state distribution can be found from the 
matrix equation Ap" = 0, and for a time-dependent part p we have a rate equation in the 
form {d/dt)p = — Ap. Using the unitary operator, U, the matrix A can be transformed to the 
diagonal form A' = U^^MA with 7^ as the diagonal elements. This transformation should 
be accompanied by the transformation of the vector p as p = Up' . Then, the vector p'{t) 
obeys the diagonal equation with a simple solution for its p-component: p'^{t) = e^'^'t^* p'^{0) . 
Correspondingly, the time evolution of the probability vector p(t) from its initial value p(0) 
is described by the formula 

{p{t))=p'-S{t)p' + S{t)p{0), (A34) 
23 



exp 



- Pa 
T 



+ 1 



where S(t) = UZ(t)U~^ , and Z(t) is the diagonal matrix with the elements Z^y{t) = 
^iivG~'^'^*. It should be noted that 5(0) = /, and S{oo) = 0, where / is the 64x64 unit 
matrix. 



7. Proton current 

The time-dependent populations, (jia) and {Np), of all redox and protonable sites in the 
model are expressed in terms of the evolving probability distribution {p{t)). We recall that 
the index a labels the redox sites a = A (Cua), L (heme a), R (heme 03), and B (Cub). The 
index f3 labels the protonable sites (3 = D,X,C, and Z. Finally, the index a labels the two 
sides of the membrane a = N, P. With the density matrix probability distributions, (p^(t)), 
over the states, of the system we can also find the proton fiows from the N-side and 
P-side of the membrane into the system, Jo- = {d/dt){N„), where {N„) is the total number 
of protons on the cx-side of the membrane: {N„) = J2q{^q(^)- Using techniques developed in 
quantum transport theory 29|, |46|], we obtain the formulas for the proton currents In and 
/p: 

In = \bD.^u\^{[l - Fn(u;^^)](p^) - Fn(u;^^)(p;,)}, 

Ip = TpJ2 \bc;>.u\^{[^ - Fp{uj,^)]{p,) - Fp{iu,^){p^)}. (A35) 

Note that these currents depend on the time-dependent probability distribution (p(t)) and, 
accordingly, they also vary with time. The total number of protons, (N^), transferred to the 
cr-side of the membrane {a = P,N) is calculated as the integral of the corresponding current: 

{N,{t))= tdhUti)- (A36) 



^0 

8. Proton-motive force 

The proton-motive force across the membrane can be defined as a difference of elec- 
trochemical potentials pp and pn involved in the Fermi distributions (IA3ip of the proton 



reservoirs: Ap = pp — pN. This gradient includes the transmembrane concentration differ- 
ence (ApH) and the transmembrane voltage V: 

Ap = V- 2.3 (RT/F) X ApH. (A37) 
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Here R and F are the gas and Faraday constant, respectively, and T is the temperature (in 
degrees Kelvin, = 1) [l, 2]- Both energy parameters, A/i and V , are measured in meV. 
At the standard conditions (T = 298 K), the concentration gradient contributes about 60 
meV per Apif-unit. This results in the transmembrane voltage V ~ 150 meV, provided 
that the total proton-motive force, A/i, is about 210 meV, and ApH = — 1 [1|]. 

The transmembrane voltage, > 0, elevates the energies of protonable sites adjacent 
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The 



to the P-side and lowers the energies of the proton sites located near the N-side 
electron sites are simultaneously experiencing the opposite effect, for the same V . As a result 
the electron energy levels, Sa-, and the proton energies, e^j, involved in the Hamiltonian Hq 

-(0). 



(lAip are shifted from their initial values, £a'* and e?^- 



ep = ef + V{xplW), (A38) 

where W is the membrane width. The positions of the redox and protonable sites, Xa and 
Xfi, are counted here from the middle of the membrane with the x-axis directed toward the P- 
side: xa ~ W/2, xl ^ xr^ xb ^ W/6, xd ~ 0.1 W, xx ~ 0.3 W, xc = 0.5 W, xz = W/6 
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FIG. 1: (Color online) Schematic diagram of cytochrome c oxidase. A single electron enters the 
enzyme at the site A and travels subsequently to sites L, R and, finally, to site B. Protons, taken 
at the N-side of the membrane, move to site which can donate protons both to the pump site 
X and to the catalytic site Z. The pre-pumped proton is transferred from the site X to the P-side 
of the membrane via site C . 
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FIG. 2: (Color online) Kinetics of electron and proton transfers in cytochrome c oxidase for 
Hp = 105 meV, ^-^ = —105 meV, T = 298 K, and V = 150 meV. The time axis (in ms) is 
shown in a logarithmic scale starting at t = 0.1 fxs. The process begins at t = 0, when a single 
electron populates the site A, and a single proton is located on the site D. (a) Time dependence 
of the electron populations ua (blue dotted line), (green dashed line), tir (red dash-dotted 
line), and ub (black continuous line), (b) The proton populations, N^) (blue dotted line), Nx (red 
dash-dotted line), Nq (green dashed line), and Nz (black continuous line), versus time, (c) The 
number of pumped protons, A'^p, as a function of time. The first phase of the process, where the 
electron moves from site A to site L, corresponds to the maximum of the population ni at the 
moment t ~ 10 fis. In the second phase, both the electron population of the i?-site, n^, and the 
proton population of the X-site, Nx, peak at t ~ 150 fis. The third phase of the process is marked 
by the significant population of the proton site Z and the electron site i? at t ~ 1 ms. In this phase 
the site X is depopulated, and the pre-pumped proton is partially transferred to the P-side of the 
membrane, A'^p 0.8. In the final phase, at t ~ 3 ms, the electron site B and the catalytic proton 
site Z are occupied, ub — 1, Nz =1, and about one proton (A'p ~ 0.95) is translocated to the 
P-side of the membrane. 
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FIG. 3: (Color online) The total energy of the system, Stot (black continuous line), and the average 
energy of the electron, E^i, as functions of time (in ms, logarithmic scale). The electron energy is 
varied in the range from E^i = —277 meV at t ~ 0, to the value Ee\ = —828 meV at t = t_B = 5 ms. 
The whole electron-proton system dissipates less energy, Ai?tot — 270 meV, than its electron 
component, which loses about 550 meV during the pumping process, indicating the energy transfer 
to the proton subsystem. 
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FIG. 4: (Color online) The efficiency of the pump, r] = Np^ts), where tB = 5 ms, as a function of 
the transmembrane voltage V at the temperatures T = 150 K (blue dashed line); T = 298 K (green 
continuous line), and T = 450 K (red dash-dotted line). For the physiological range of transmem- 
brane voltages, 150 meV < V < 200 meV, the pump operates with a maximum efficiency at room 
temperatures (T = 298 K). At high temperatures and high enough transmembrane voltages, the 
efficiency r] takes negative values, suggesting that at these conditions the protons flow back, from 
the positive to the negative side of the membrane. 
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